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We study polarized Semi-Inclusive Deep Inelastic Scattering (SIDIS) processes, £(Se) + p(S) — > 
£' hX, within the QCD parton model and a factorization scheme, taking into account all transverse 
motions, of partons inside the initial proton and of hadrons inside the fragmenting partons. We 
use the helicity formalism. The elementary interactions are computed at LO with non collinear 
exact kinematics, which introduces phases in the expressions of their helicity amplitudes. Several 
Transverse Momentum Dependent (TMD) distribution and fragmentation functions appear and 
contribute to the cross sections and to spin asymmetries. Our results agree with those obtained 
with different formalisms, showing the consistency of our approach. The full expression for single 
and double spin asymmetries As e s is derived. Simplified, explicit analytical expressions, convenient 
for phenomenological studies, are obtained assuming a factorized Gaussian dependence on intrinsic 
momenta for the TMDs. 

PACS numbers: 13.88. +e, 13.60.-r, 13.85.Ni 

I. INTRODUCTION 

Experiments with inclusive Deep Inelastic Scattering (DIS) processes, £N —> £' X, have been performed for 
decades and have been interpreted as the most common way to investigate the internal structure of protons 
and neutrons. At large energy and momentum transfer the leptons interact with the nucleon constituents; by 
detecting the angle and the energy of the scattered lepton one obtains information on the partonic content of 
the nucleons. This information is encoded in the Parton Distribution Functions (PDFs) which give the number 
density of partons moving collinearly with the nucleon and carrying a fraction x of its momentum at a certain 
value of the squared momentum transfer Q 2 . The prediction of the Q 2 dependence of the PDFs has been one of 
the great successes of pQCD. Although successful, such an approach only offers information on the longitudinal 
degrees of freedom of quarks and gluons, giving no information on the transverse motion, which is integrated 
over. This transverse motion - transverse with respect to the parent nucleon direction - is related to intrinsic 
properties of the partons, like orbital motion, and reveals new aspects of the nucleon structure. 

In the last years, driven by unexpected spin effects and azimuthal dependences, the study of the intrinsic 
motion of partons has made enormous progress; indeed, a new phase in the exploration of the proton and 
neutron composition has begun. The leading role in such an effort is played by Semi-Inclusive Deep Inelastic 
Scattering (SIDIS) processes, IN — > £' h X, in which, in addition to the scattered lepton, also a final hadron 
is detected; this hadron is generated in the fragmentation of the scattered quark (or gluon) - the so-called 
current fragmentation region - and, as such, yields some new information on the parton primordial motion. 
This new information is encoded in the so-called Transverse Momentum Dependent partonic distribution and 
fragmentation functions (TMD-PDFs and TMD-FFs, or, shortly, TMDs), f a/p (x,k±) and D h/a (z,pj_). The 
TMD-PDFs give the number density of quarks (a = q) or gluons (a = g) with light-cone momentum fraction x 
and transverse momentum k± inside a fast moving proton; the TMD-FFs give the number density of hadrons h 
resulting in the fragmentation of parton a, with a light-cone momentum fraction z and a transverse momentum 
Pj^, relative to the original parton motion. At leading- twist, taking into account the parton and the nucleon 
spins, there are eight independent TMD-PDFs [2 [5]; if the final hadron is unpolarized or spinless, say a pion, 
there are two TMD-FFs. All these quantities combine into physical observables and by gathering information 
about them one accesses the momentum distribution of partons inside the nucleons. 

The theoretical framework used to analyze the experimental data is the QCD factorization scheme, according 
to which the SIDIS cross section is written as a convolution of TMDs and elementary interactions: 

fotr+fhx = ^ / 9/p ( Xj fe±; q2) g, dfj ^ g JJ h/q (z,p x ; Q 2 ) . (1) 
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In the 7* — p cm. frame, see Fig. 1, the measured transverse momentum, Pt, of the final hadron is generated 
by the transverse momentum of the quark in the target proton, k± , and of the final hadron with respect to the 
fragmenting quark, p^. At order fcj_/Q it is simply given by 

P T =zk x +p x . (2) 

There is a general consensus [3H7] that such a scheme holds in the kinematical region defined by 

P T ~ Aqcd « Q ■ (3) 

The presence of the two scales, small Pt and large Q, allows to identify the contribution from the unintegrated 
partonic distribution (Pt ~ k±), while remaining in the region of validity of the QCD parton model. At larger 
values of Pt other mechanisms, like quark-gluon correlations and higher order pQCD contributions become 
important [THS]. A similar situation [4l l6l IT0HT6] holds for Drell-Yan processes, AB — > £ + £~X, where the two 
scales are the small transverse momentum, qT, and the large invariant mass, M, of the dilepton pair. 

Let us elaborate now on Eq. 0. We consider the SIDIS cross section at the leading a em order - i.e. one- 
photon exchange - and in the "standard" |17) kinematical configuration of Fig. 1, which defines the azimuthal 
angles <j>h and </>s in the 7* — p cm. frame. The most general dependence on these angles has been discussed 
in several seminal papers [Tl [T8H2"U] . both in a model independent scheme and in the parton model. According 
to the usual derivation, the polarization states of the virtual photon, as emitted by the lepton in a certain 
direction, contains azimuthal dependences |18l 119] ; within the parton model, the virtual photon scatters off a 
quark - which subsequently fragments into the final hadron - and each term of the azimuthal dependences can 
be written as a convolution of distribution and fragmentation functions (TJ ITM2"2"] . 

We re-derive here the same general expression of the cross section, and its parton model content, by assuming 
from the beginning the validity of the TMD factorization ([lj ; we use the helicity basis to compute the elementary 
interaction and to introduce transverse momentum dependent distribution and fragmentation functions. In such 
an approach the full azimuthal dependence is simply generated by the properties of the helicity spinors and 
amplitudes. Our final results coincide with the existing ones, showing the full equivalence of the two procedures. 
Our formalism is based on a physical and intuitive picture, which somehow factorizes the physical process in 
different steps: the "emission" of a parton by the interacting hadron (p — > q + X), the interaction of the parton 
with the lepton (Iq — > Iq), and the "emission" of the final hadron by the scattered quark (q — s- h + X); each 
step is described by the corresponding helicity amplitudes. For SIDIS processes this factorization has been 
formally proven and expressed in terms of TMDs, Eq. (JlJ. Such a procedure can naturally be extended to other 
processes, and indeed this has been done for the large Pt production of a single particle in inclusive hadronic 
interactions, AB — > C X [2J. The point, however, is that, despite the natural simplicity of the approach, the 
TMD factorization has not been proven for processes with a single large scale, like A5->CX. Due to this, the 
study of dijet production at large Pt in hadronic processes was proposed [23l426j . where the second small scale 
is the total qT of the two jets, which is of the order of the intrinsic partonic momentum k±. This procedure 
leads to a modified TMD factorization approach, with the inclusion in the elementary processes of gauge link 
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color factors [27-30 . However, some doubts on the validity of such a factorization scheme have been recently 
cast [3T] . A possible experimental test of the TMD factorization for processes with only one large scale has been 
proposed in Ref. [32] • We limit our discussion in this paper to SIDIS processes, in the kinematical region g) for 
which TMD factorization holds, and obtain the most general expression for the polarized cross section, with our 
hclicity formalism. A similar study can be done, with the same validity, for Droll- Yan processes [T2J [23 |33] . We 
introduce only leading-twist TMDs and take into account exact kinematics, often simplifying results by only 
keeping terms up to 0(k±/Q). 

The paper is organized as follows. In Section II we present our formalism and compute the polarized SIDIS 
cross section. In Section III we give the explicit general expressions of all independent single and double spin 
asymmetries, in terms of the TMDs. In Section IV we give explicit analytical formulae for the spin and azimuthal 
asymmetries, assuming a factorized Gaussian dependence of the TMDs on k± and p±. In Section V we draw 
our conclusions. Useful results are derived and collected in Appendices A-E. 



II. CROSS SECTIONS IN POLARIZED SIDIS 



According to Refs. [34] and [2] the full differential cross section for the polarized SIDIS process, £(St) +p(S) — > 
£' hX, can be written, within TMD factorization, as 

d(J i{s e )+ P (s)^i'hx _ 1 i r z_ 

dx B dQ 2 dz h d 2 P T d0 s ~ ^gf^ } 16 7r(x B s) 2 7 ± z h 

x ^l^C W . (4) 

where we adopt the kinematical configuration of Fig. 1, and, as usual: 

■s = (i + P ) 2 Q 2 = -q 2 = -(£-£') 2 x B = ~ = — ■ ( 5 ) 

2p ■ q p-q 

The variables x, z and p ± which appear under integration in Eq. Q are related to the final observed variables 
x B , Zh and Pt and to the integration variable k±. The exact relations can be found in Ref. |34| : at 0(k±/Q) 
one simply has 

x = x B z = z h p ± = P T z h kj_ . (6) 
J includes some non-planar kinematical factors [34) : 



x B ( \ xi k 
x 



2 1,2 v- 1 



where the last relation holds at 0(kj_/Q). At this order Eq. Q can be written as: 

da e(s e )+ P (S)^e'hx 



dx B dQ 2 dzh d 2 Px dtps 2 



1 {A} 



where we have explicitly shown the integration over p ± for clarity and further use. In Eqs. Q and ([8]) the sums 
are performed over all quark flavors (q — u, u, d, d, s, s) and all quark, lepton and hadron helicity indices; p^ S y 
is the initial lepton helicity density matrix, which describes the spin state of the lepton beam; for unpolarized 
lcptons one simply has p e x A , = ^S\ e \' e - It might be helpful, and useful for physical interpretations, to recall 
that, in general, for a spin 1/2 Dirac particle one has: 

Pxy ~ 2 I P x + iP v 1-P, (9J 



where Pj = P x ,P y ,P z are the components of the particle polarization vector in its helicity frame (throughout 
the paper we follow the definitions and conventions for helicity states of Ref. [35 ) . 
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Let us discuss in detail the different "factors" in Eq. Q : they represent the distribution of polarized partons 
(only quarks at LO) inside the proton, their interaction with the lepton and the fragmentation of the (polarized) 
final quark into the observed unpolarized hadron h. We follow, and adapt to the case of SIDIS, the discussion 
of Ref. P]. We describe the three stages of the process - quark emission, interaction and fragmentation - within 
the helicity formalism, which allows us to introduce in a natural way, at each step, several phases; these, when 
combined into the expression for the physical cross section Q give its full azimuthal dependence, in agreement 
with results in the literature derived in a more formal and somewhat less intuitive way |22j . 



A. TMD partonic distribution functions 

P\^ P }f fq i /p,s{ x T^ t '±) counts the number of polarized quarks inside a polarized proton; it is the polarized 
distribution function of the initial quark qi with light-cone momentum fraction x and intrinsic transverse mo- 
mentum kj_, inside the target proton p in a pure spin state S. Using Eq. ([9| and parity invariance one can see 
that there are eight independent distribution functions, which can be defined as: 

P?f q/p ,sA*M = /« /s >,feL) - fl Sj/ST (^k ± ) = A/« /ST (z,fc x ) (10) 
Pj f q / P ,s L (x, k ± ) = f% /SL (x, k ± ) - f\ j/SL (x, k ± ) = A/« /Si (x, fex) (11) 
fg/ P ,s T {x, k±) = f q/p (x, k±) + - Af q/ST (x, k±) , (12) 



with 



We define, for further use, 



Afq/s T (x, k±) = f q /s T (x, k±) - fq/s T (x, k±) . (13) 



\ [fs y /s T (x,k±) - f Sy /s T (x,k ± )} ee A f Sy/ST (x,k ± ). (14) 



In Eqs. (10) and (111, j = x,y,z are the coordinate-axes in the quark helicity frame and Sl,t are respectively 
the longitudinal and transverse components of the proton polarization vector, with respect to its direction of 
motion. 

Different notations can be found in the literature for these functions, in particular those introduced by the 
Amsterdam group 1, 36, 37 , which are largely adopted. The relationships between the two sets can be found 



in Ref. [2], and will be repeated for convenience in Eqs. ( 22 )— ( 25 ) 



According to the physical interpretation of the factorization scheme, as outlined above, these quantities can 
be introduced by making use of the helicity amplitudes F\ ,a x ;A , which describe the soft process p — > q + X. 
Since the partonic distribution is usually regarded, at LO, as the inclusive cross section for this process, the 
helicity density matrix of a quark q inside the proton p with spin 5* can be written as 

q/p,S f / , % _ v,s 



V.A.v q '" X ' K 



having defined 



'tv Z ^K^kH M ': (16) 



pi p 



A , A > 



where the stands for a spin sum and phase-space integration over all the undetected remnants of the 

Jx,x x 

proton, considered as a system X, and the .F's are the helicity distribution amplitudes for the p — ¥ q + X process. 



Eq. ( 15 ) relates, via the unknown distribution amplitudes, the helicity density matrix of the parton q, 



q/p ,s _1 1 + PS PS - iP$ \ _ I ( 1 + P q z P q T i n 7 ^ 

q^'q 2 V P% + iP§ 1 - -PI J 2 v P| e iv ><> 1 - PS )' 



5 



to the helicity density matrix of the polarized parent proton, 



1 + Sz Sx — iSy 
Sx + iSy 1 — Sz 



1 + S L S T e-^ s 
S T e lvs 1 - S L 



(18) 



In the above equations S — (Sx,Sy,Sz) — {St cos (ps, St simps, Sl) is the proton polarization vector and ips 
its azimuthal angle, defined in the helicity reference frame of the proton p. Similarly, P q = (P£, P q , P q ) = 
(Pj, cos (p s , P T sin tp s , P!*) is the quark polarization vector defined in the quark helicity frame and ip Sg its 
azimuthal angle. For the kinematical configuration of Fig. 1, one has (ps — 2tt — (f>s (see Appendix |B| , so that: 



p.S 



S 



1 + Sl &t < 

its 1 _ s 



(19) 



Notice that, in general, we denote by ip angles defined in the proton or quark helicity frames, while the symbol 
<t> is used for the corresponding angles measured in the 7* — p cm. frame. 

The distribution amplitudes F depend on the parton light-cone momentum fraction x and on its intrinsic 
transverse momentum kx , with modulus kx and azimuthal angle <px , in a precise way [21 135] , which, again 
referred to the kinematical configuration of Fig. 1, reads: 



Fx q ,x x -x p (x,k ± ) 



so that 



$(x, kx) = F^(x, k ± ) exp[i(\' p - A p )0i 



(20) 



(21) 



F x q \?{x,kx) has the same definition as F x q 'y(x, fej_), Eq. iWh with F replaced by F, and does not depend 
p'p p : p 1 1 

on phases anymore. Notice that we have chosen, throughout the paper, to denote with a hat all soft quantities 
which depend on both the modulus and the phase of the kx and intrinsic momentum vectors, while we drop 
the hat for quantities which only depend on the modulus of these vectors and not on their phases. 

Eqs. (15 1, (17), (19) and (21), together with parity properties and the arguments collected in Appendix B, 
allow to extract the explicit phase dependence of the eight polarized distribution functions (10)-(12), with the 
result (more details can be found in Ref. [2]): 



f q/PtS (x,k ± ) = F++{x,k 1 _)+F++(x,k ± )-2S T lmF++(x,k ± )sm(cf >s -<l> 1 _} 



(22) 



f q / p (x,k ± ) 
fi(x,k±) + 



■ St Af q/ST (x, k±) sin((f>s - <j)x) 



k± 

St -jt f 1T {x,kx) sin(0 s 



PI ? q / P ,s( x > k ±-) = Sl [Ft+( x > ~ F++{x, k x )] + 2 S T ReF++(x, k±) cos(0 s - <t>±) 

+ i 



S L AJI /Sl (x, kx) + S T A/« /Sr (x, k ± ) cos(^ s - <f> x ) 



Sl 9il{x, kx) + S T ^ 9rr(x, kx) cos(</> s - 4>x) 

P q x f q / P ,s(x,kx) = -2S L RcF++{x,kx)-S T [F+I(x,kx) + F+ + (x,kx)] cos(0 s - (j>x) 
= -S L Afl /SL (x,kx) - S T Afl /ST (x,kx) cos(^ s - <t>x) 



(23) 



(24) 



kx 

-Sl h 1L (x,kx) - S T 



hi(x, kx) + 



k 2 

2M 2 hlT 



(x, kx) 



cos 



-<f>±) 



P« f qM s(x,k x ) = 2lmF+-(x,k ± ) + S T [F+Z(x,k ± )-F-+(x,k ± )] sin(<fe - 4>x) 
= -Af« v/p (x,kx) + S T A-fl /ST (x,kx) sin(0 s -^) 



(25) 



kx 

M 



h x {x, kx) + St 



hi(x,kx) 



2M 2 



hi T (x, kx) 



sin 



As already stated, 4>s & n d 4>x are respectively the azimuthal angle of the proton polarization vector S and of 
the quark intrinsic momentum k± measured in the 7* — p cm. frame of Fig. 1. Also the quark polarization 
vector components P? (i — x, y, z) refer to the helicity frame of the quark, as reached from the 7* — p frame: 



G 



this explains the sign differences between Eqs. ([22} |24}|25) and Eqs. (B12, B14-B15) of Ref. (in the latter 
case the polarized proton was moving along Z cm rather than — Z cm . Further comments are given in Appendix 
B). Notice that, while P« f q/p ^ 0, one has P« f q/p = PJ f q/p = 0. 

The above equations, which will be soon used, deserve some further explanation. In each equation the first line 

expresses the partonic distributions in terms of the F^ q '^ (x, fcjj's and shows their exact phase dependence. The 
second line gives the same quantities using our notations for the TMD-PDFs. According to our "hat convention" , 
quantities like A/'.^ s (x, fcjj do not depend on phases anymore, as such dependence has been explicitly extracted 
out; comparing with Eqs. ( 10 )— ( 12 ) one has (always referred to the variables and kinematical configuration of 
Fig. 1): 



Af q /s T (x,kx) = -Af g /s T (x, kx) siring - 4> x ) 
Afl /SL (x,k ± ) = -Afl /SL {x,k ± ) 

Af^ /ST (x,kx) = -Af^ /ST (x,k±_) cos(cf)s - <f>±) 



Af Sy/SL (x,k ± ) = 
Af« y/ST (x,k ± ) = 



-Afl / „(x > k ± ) = -Afl,(x,k ± ) 



f'l 



_J p \K,^> "I " J Sy /S T 

-A/ s 9 H/p (x,^) + A-/ s ' y/Sr (x,fc ± ) 



■ilv 



x) sin 



A/! 



/s L 



(x,k ± ) = Af^ /SL (x,k±) 



(26) 
(27) 
(28) 
(29) 
(30) 

(31) 
(32) 



A fl/s T ( x ' fc ->-) = AJ^ /St (x, k ± ) cos(c/) S - <j> ± ) . 

According to our choice the Af^ , Srr Sl (x, k±) introduced here are the same as in Ref. [2]. 

The last line of Eqs. ( 22 1— ( 25 ) gives the connection with the Amsterdam group notations; M is taken as 
the proton mass. These last relationships hold at leading twist; notice also that, when comparing with the 
results of the Amsterdam group, one should take into account other differences in conventions and notations. 
In particular: 



(.Pt ) Amsterdam 
( — z kx) Amsterdam 



k ± 

Px = {Pt 
Pt 

Finally, we recall some other notations widely used in the literature 



(h) 



Amsterda 



2k ± 



A N f q/p ,(x,k ± ) = Af q/ST (x,k ± )=41mF+±(x,k ± ) = -^/^(i.fci) 



A N f q y p (x, k x ) = Af q s , (x, kx) = -2ImP++ (x, k±) 



M 



(33) 
(34) 

(35) 



(36) 
(37) 



1 

2 

l r 

2 



a/? 



,(x,kx) + A-f q Sy/ST (x,kx) = F+_(x,kx) = h 1T (x,kx) + 



^hi T (x,kx) = h 1 (x,k ± ) (38) 



(x, kx) - A-f« y/ST (x, kx) = F++(x, k ± ) 



k 2 

— ^rhi~r r (x, k±) 
2M 2 1TV ' 1 



(39) 



Axq{x) — hi(x) — Id kx hi(x, kx) = d kx 



h 1T (x, kx) + 



2A'p hlT 



(x,kx) 



(40) 



Eqs. (36), (37) and (40) refer, respectively, to the Sivers, the Boer- Mulders and the transversity distributions. 



B. TMD fragmentation functions 

The quantity D x h ' w (z,Pj_) describes the hadronization of the quark qr into the observed final hadron h, 
which carries, with respect to the fragmenting quark, the light-cone momentum fraction z and the intrinsic 
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transverse momentum p x . Similarly to the distribution functions, also D^ h 'y h (z, p_|_) can be written as the 
product of fragmentation amplitudes for the q — > h + X process: 



(41) 



where the J stands for a spin sum and phase space integration over all undetected particles, considered 

as a system X. The usual unpolarized fragmentation function D h / q {z), i.e. the number density of hadrons h 
resulting from the fragmentation of an unpolarized parton q and carrying a light-cone momentum fraction z, is 
given by 



D 



h/q 



(42) 



We consider only the cases in which the final particle is either spinless (A h = 0) or its polarization is not 
observed, 



(43) 



In such a case, parity invariance reduces to two the number of independent D*\\' (z,pj_). These, in general, 
may depend on the azimuthal angle of the final hadron momentum 1\ around the direction of the fragmenting 
quark q, as defined in the quark q helicity frame, which we denote by ip q (it was actually denoted as (j)q m 
Ref. 0): 



D h 4(z,p x ) = D h Jl{z, Pi _) = D h/q (z, P± ) 

D>l / l(z,p ± ) = D h + / «(z,p 1 _)e^ 

D h 4(z, Pl _) = [D'^(z,p ± )Y = -D h ^(z, Px ) e~< 



(44) 
(45) 
(46) 



In Appendix [c] it is shown how to express ip q in terms of integration and external variables (defined in the 
7* — p cm. frame), with the result, at leading order in the (k±/Q) expansion: 



COS ipg 



Pt 

p± 



cos (4> h -4>±.)- z h 



Pt 



h p T 
sin<^ = — sin(0h - cj>±) . 



(47) 
(48) 



In Eq. (44) D h / q (z,p±) is the unintegrated unpolarized fragmentation function. Other common notations used 
in the literature are: 



A N D h/q r(z,p ± ) ee -2iD h 4( Z , Pl _) = 2lmD h + ^(z, Pl _) = ^-H^(z, P± ) 



(49) 



referred to the Collins fragmentation function. Mh is the mass of the produced hadron. 



C. Elementary interaction 

The M\ t , \ q ; a £ a,. are the helicity amplitudes for the elementary process tq% —¥ £'qf, computed at LO in the 
7* — p cm. frame, taking into account the quark intrinsic motion; the amplitudes are normalized so that the 
unpolarized cross section, for a collinear collision, is given by 

d& e qi ^e' qf ii 

= ie^iS'^/^l ' (50) 



where t = —Q 2 and s = x B s. 
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Helicity conservation for massless particles requires A^ = X^, X q . — X q , = A g , which implies that there are 
only two independent non- vanishing amplitudes, explicitly computed in Appendix [Aj with the result: 



where y 



Mi = M++.++ = Ml_.__ = e q e 2 
M 2 = M+_ ;+ _ = M*_ + ._ + = e q e 2 
- and 



- A+ e +l 

y 

i-y 



y 



A + e 



i-y 



i 



4 ^ 1 z y 

y Q 
y y Q 



(51) 
(52) 



A+ = l±Wl + 4-i 



(53) 



These are exact LO results, holding at all orders in the k±/Q expansion. By truncating this expansion at first 
order in k±/Q, one obtains much simpler expressions, which will be useful later, 



Mi =M ++;++ ~2e,e 2 



1 



Ma = M 4 



2e„e' 



(l-») 



-i<t>± 



yi - y k±_ 
y Q 

v Q 



(54) 
(55) 



We can now assemble the expression of the different factors - each corresponding to a physical step - into 
Eqs. Q or pj) to obtain the SIDIS cross section in terms of the TMDs. This can be done in several ways. The 
most direct one is that of performing the helicity sums in Eq. Q taking into account Eqs. (17), ( 44 1— ( 46 ) , (49), 
d5il and d52|. It yields: 



d(J e(s e )+ P (s)^£'hx 
dx B dQ 2 dzh d 2 PTd(j)s 



-E 



2-7T ^— ' 16 7T (x B s) 2 



d z kj_ — J 

Zh 



\ j/,/p,s(a;,fe±) (|M!| 2 + |M 2 | 2 ) D h/q (z, P± ) 



+ P*P«f q/p>s (x,k ± ) [iM^-lMtf) D h/q (z, Pl _) 



(56) 



PS f q / P A x > (Re(M x M 2 *) cos $ ~ ImiMiMi) sin ^ 
P« f q/p>s (x,k x ) (MMiM^) cos ^ + Re(M 1 M 2 *)sin^) A N D h / q -r(z,p±] 



which expresses the cross section in terms of the lepton and the quark polarization vectors, the helicity am- 
plitudes of the elementary interaction and either the unpolarized or the Collins fragmentation functions. The 
intrinsic transverse momentum of the produced hadron, p^ , is related to k± and the other kinematical variables 
as shown in Eq. (28) of Ref. [34]. The exact expressions of cos ip q and sin<^ can be obtained from Eqs. (C3) 
and (C4). 



We now continue our computation, in this Section, at 0(k±/Q). From Eqs. (54), (55), (47) and (48), we 
have: 



\M 2 



4e„e 



2 C 4 



\Mi I 



Aeie 



1 + (1 - yf - 4(2 - y) y/l^y ^ cos cf> ± 

k x 



1 - (1 - j/) 2 - AyyT 



y -q cos cp± 



lm{M x M^) cos^' 1 + Re(M x M 2 *) siny£ = ^L^_ J (i _ , n 



p± y 



sm{<p h 



^— sm20j_ 
"t 



Q 



sin h - z A — sin ± 



(57) 
(58) 

(59) 
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Re(MiM 2 *) cosy?* - Im(MiM|) sin^ 



P r 4e 2 e 4 

p± y 2 



(i-y) 



2^T^j(2-y)-± 



cos(0/, + 4>±) - Z h — COS20_L 



COS - Z ft — COS 



(60) 



Inserting these results, together with Eqs. (22|-(25|, into Eq. (56), gives, at order k±/Q, the following 



expression for the SIDIS cross section in the TMD factorization scheme: 



dx B dQ 2 dzh d 2 Px d(f>s 2tt 
\f q/p [1 + (1 - yf] D h/q - Uf ^(1 - y) 



P± 



cos(<f> h + <t>±) - z h — COS 20_L 



4e 2 e 4 

y 2 



-2(2 -y)v^y 



1 P T 

f q/p coacf> ± D h/ g--Af q Sv/p — 



P-L 



(1 - y) A ft/s L (s' m ((t>h + <(>±) - Zfc ^ sin20 ± ^ 



P T 

h/qT 



coscj) h - z h cos0_l ) A N D h/q 



-2(2 - y)^T^j^L ^ a/« -/5i (sin0 h - % ^ sin0 ± J A%, t 
f [! - (1 - y) 2 ] W./Sr.Dh/, ~ 4yv^ ^ A/« /S£ oos^£» 



Q 



- [1 + (1 - yf] Af q/ST sin(0 ± - 4>s)D h/q 



+Pj [1 - (1 - yf] Afl /ST cos(0 ± - <j> s )D h/q 



-P^ 2yyT^ A/« /St ( cos S + cos(20 ± - 0s)) A*/, 



+ 



^ (1 - y) (A/ S 9 x/Sr + A-^ /Sr )(sin(0, l + S ) - z„ ^ sin(0 ± + S )) A N D h/q , 



2p 

+ (1 ~ V) (A/ ^/ s - " A_ - f l/s r )( sin (^ + 2 ^ ~ - * h |^sm(30 ± - <i> s ))A N D h/q , 
(2 - y) v /T3^|L ( A/^ /ST + A-^ /ST )(sin(^ - X + 0s) -z h ^- sin0 s ) A" D h/q 



Pi 
Pt 



Q 



--^ (2 - y )^l^-±(Afl /ST - A-/« s/Sa .)(Bin(^ + 0x - S ) - z fc j± sin(20 ± - S )J A^ t/ , t 



+(2 - y)v / T 3 y ^ A/, /St (sin 05 - sin(20^ - </>s))D h/q J . 



(61) 



The first three terms of Eq. (61 ) correspond to the contribution of the unpolarized proton to the SIDIS cross 
section; they contain either the unpolarized or the Boer-Mulders distribution functions. The following three 
terms correspond to the longitudinally-polarized proton contributions; they depend either on the helicity dis- 
tribution A/* [= Aq = gi] or on the A/' , Sl [— (k±/M) hj^ L ] transverse momentum dependent distribution. 
Finally, the last eight terms correspond to the transversely-polarized proton contributions; they may originate 
from the Sivers function, from A/' / St [— (k±/M) gyp], and from the transversity distribution functions, related 
to the combinations (A/' ,„ ± A - /' ,„ ) as shown in Eqs. (38) and (39). The partonic distributions couple 
either to the unpolarized or to the Collins fragmentation functions, depending on whether they are, respectively, 
chiral even or odd. 

Notice that we have intentionally grouped all terms according to their phases, so that this expression can 
be easily compared with the analogous formulae of Ref. 22J, which have the same structure. To make the 
comparison fully explicit, apart from converting our notation to the Amsterdam group notation, we need to 
extract from the integration over the intrinsic transverse momentum k± the dependence on the azimuthal angles 
4>h and 4>s- On the basis of a simple tensorial analysis, which is described in detail in Appendices [D] and [E{ we 
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can recover Eqs. (4.2)-(4.19) of Ref. [52], without formulat 
(p_i_) dependence of the distribution (fragmentation) functions. 

In analogy with the Amsterdam notation, Ref. [22], we define the convolution on transverse momenta in the 
following way 

C[wfD]=J24 

9 

Notice that this definition differs from Eq. (41) of Ref. [22 by a factor x B and 

i of Eq. (61 1 can in fact be written as 



can recover Eqs. (4.2)-(4.19) of Ref. [22], without formulating any particular assumption on the x (z) and k± 
(px) dependence of the distribution (fragmentation) functions. 

In analogy with the Amsterdam notation, Ref. 22J, we define the 
following way 

J d 2 k ± d 2 p ± 5^(P T -z h k ± -p ± )w(k ± ,P T )f{x B ,k ± )D(z h , P± ) . (62) 

for the definition of the parton 



iNoticc tnat tins acmution dim 
momenta, see Eqs. (33)-(35). 

The convolutions on intrinsic transverse 



Fuu = 



momenta in the single terms 
Y4 f d 2 kxf g/p D h/g = C[f 1 D 1 ] 



cos l<pb t vu 



-Ye 2 a [ d 2 k ± Af , ^ 
Z_, i J Js «/P2p 1 _ 



cos(0/j + 0_l) - z h — cos 2(j)±_ A N D h/q1 - 



— cos 2cj)h C 



(P T ■ k±) - 2z h (P T ■ kx) 2 + z h k\ . , 
7. A,f, U ft i 



COS (phtjjjj — 



-2J24 l^^jcos^ 
Pt 



z h M h M 

Jq/p D h/q 



= COS 0k 



2 s 



c 



2 P± 
{Pt ■ k±)ft Di 



COS <ph - Z h — COS0j_ 



Af Sy/p A N D h/q ,} 



k 2 ± (p t -z h P T - kx 



z h M h M 



^ Hi 



sm2<p h F UL Vh 



Ye 2 [ d 2 kxP^Af q /q 



2p ± 

■kx)-2z h (P T -kx) 2 + z h ki ± j 
z h M h M 1L 1 



kx \ 
x)- z h — sin 2<j) x j A N D h/q r 



(Pt ■ kx 



sin 0^ i' (7L 



sin 2<j) h C 

- 2 Y. e l J d2fe± ^27T A/ ^/^ ( sin ^-^^ sin ^J AJV£ W 



= sin^ 



k\ 



P T - z h (P T ■ kx)) 
: ^T^f 1 



z h M h M 



T sin 0^ 



sm <p h F LU 



Fll 



cos (j> h F LL 



.cos 4> h 



= at leading twist 

= Y e ll <? k ± A /^/5, ^v, = c iaiL Dx] 

q J 

d2k± ~$ A ft*/St. C0S( f>xD h/q 
(P T ■ kx)gi L D 1 

t)X - 4>s)D h /q 



= COS(/) h 



-De 



Mth-^F^-^ = \Y e l Jd 2 kxAf q/ST sin( 



= sin(0/, -<j>s)C 



M JlT 1 



cos^-^F^-^ =Y,< I d 2 kxAfl /ST cos(cf>x-<P s )D h/q 

9 J 



(63) 



(64) 



(65) 



(66) 



(67) 

(68) 
(69) 



(70) 



(71) 
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.COS (fis 



COS Cps b LT 



= cos(<f> h - (/> s ) C <kPT ] ^ 1 ^ 9lT D l 
= - I d 2 k^Afl /ST cos<j> s D h/q 



= cos0 s f-Ij C 



(72) 



(73) 



cos(2^-0 5 )i^ T s(2 ^ s) 



= cos 



fci - 2(P T • fc ± ) 2 ) 



1 



{24> h ~4>s)-^C 



M 



9lT D l 



(74) 



k 



( A fl/s T + A 'fl/s T )(M4>h + <Ps) ~ z hpr 

' (p T - z h k±{P T • k±j) 

z h M h 11 



sin(0j_ 



-4> s j)A N D h/qt 



= sin((f> h + <j) S ) C 



sm{6<ph - q>s)r UT 



(75) 



= W [d 2 k 



(Af q sJST - A-/^ /Sr )(sin(^ + 20 ± - S ) - z h ^- sin(30 ± - s ))A^Z» h/tft 



= sin(30;, - S ) C 



■ 2 ± {-P T + 2P T (P T -k ± ) 2 -z h k ± [4(P T -fe ± ) 3 + 3(P T -fc ± )] } 



2z h M h M 2 



foyp Hi 



(76) 



Pt k± 

~p± Q 



<(Afl /ST + A-/ s 9 y/ST )(sin(^ - ± + S ) - 2^ sin0 s ) & N D h/qt 



1 f k 

2^2 e l J d2fe -L-^ A /g/ST sin 0s ^/-j 



= sin 0s 



sin(20 h - 0s)P^t 



• 2 d 2 kl £ZL*± 



(Pr-^- Zh kl) ^ 

z h M h 1 2M JlT 

Pt k± 



(77) 



x(A.f^ /Sr - A-/ s 9 y/5r )(sin(0,. + X - S ) - z,,^- sin(20 ± - S )) A^.^t 
~ \Y.< J d 2 k 1 k -^Af q/ST sin{2^- < Ps)D h/q 

fci((P r -fc ± ) + z /l fci(l-2(P T .fe ± ) 2 )) ± 

hlT ^1 



= sin(20 h -0 s ) (--) C 



Mh M 2 



Zh M h 



(2{P T -k 1 _f-k 2 1 ) 
- ~ /it -Di 



(78) 

These "Fs e S structure functions" are the same as those defined in Ref. [52], apart from an overall factor x B 
which appears in the latter. In the comparison one should consider only leading twist TMDs and remember the 
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Eqs. ( 33 1— ( 35 ) . Using the above F's in Eq. (61) one obtains the full expression 



different notations of Ref. [ 
of the SIDIS polarized cross section, valid with leading twist TMDs and at kinematical order k±/Q: 



da e(s e )+ P (s)^e'hx 
dx B dQ 2 dzh d 2 Pr dcfis 



2 a 2 



1 + (1 2 V? F vv + (2 - y)^~y cos <f> h + (1 - y) cos 2<t> h F™ 2 ^ 



+ S L 



S L P : 



(1 - y) sin 20, F$?*» + (2 - y)y/l^y sin0, F^ 

i-(i-y) 2 



f L l + yy/i-y cos <f) h F^ ' 



+ S T 



l + {l-yf 



, \ J? sin(0 h -<^ s ) 



+ (l-y) (sin(0, + S ) F ^+^> + sin(3 ^ _ ^ ^ 
+ {2-y) Jl^y ( sin S F^ s + sin(20 h - <j> s ) 
l-(l-2/) 2 



,sm(3<j>h— <t>s) 
UT 



sm{2(f> h —<f>s) 



COS(<Ph 



L \ TpC°s(4>h-<t>s) 



+ 



y V 1 - V ( cos F^t 05 +cos(20, i - (f) S ) 



f: 



cos(2<f>h — <l>s) 



LT 



(79) 



This expression agrees with Eq. (2.7) of Ref. [22], bearing in mind Eqs. (2.8-2.13) and that, at leading twist, 

*uu,l - f LU - U. 

In obtaining the general cross section structure of Eq. ( 79 1 we started from the TMD factorization, Eq. Q ; 
then we have simply exploited the properties of the helicity amplitudes, which essentially originate from the 
phase dependence of the Dirac spinors and their non collinear kinematics. Each step of the factorization scheme 
contributes some phases, including the elementary interactions. 

Some of the final azimuthal dependences have a clear and direct physical interpretation. For example, 
the phase of F^ T ^ h ~^ s \ Eq. (71), originates from the phase dependence of the Af q / ST (x,k±) distribution, 
Eq. ( 26 ) . This is the Sivers effect [551 [3"9~| , which relates the number of unpolarized quarks with intrinsic 



momentum fcj_ to the spin of the proton; such an effect, due to parity invariance, can only be of the form 



S ■ (px k±) = St sin(0j_ — 0s). Similarly, the phase in the first term of F^ , Eq. (65 ), being associated with 



unpolarized distribution and fragmentation functions, can only come from the kx dependence of the elementary 
interaction, the so-called Cahn effect 1531. 



III. SINGLE AND DOUBLE SPIN ASYMMETRIES IN SIDIS 



From the expression of the SIDIS polarized cross section we can now compute all spin asymmetries which 
have been, or can be, measured. We can restart from Eq. (56), inserting into it the expressions of the polarized 
quark distributions, as given in Eqs. ( 22 f— ( 32 ): 



da l{S t )+ P (S)^i'hX 



Ly L 



dx B dQ 2 dzh d 2 PT d(j)s 2n ' 16 ir (x B s) 2 



d z kj_ — J 

Zh 



1 

X - 

2 



f q/p (x,k x ) + -S T Af q/ST (x,k ± ) ) [\M L \" + \M 2 \ 2 ) D h/q (z, Px ) 



P l z (S L Afl /SL (x,k ± ) +S T Afl /ST (x,k x )) (|Mx| 2 - |M 2 | 2 ) D h/q (z, P± ) 
Af q sy/p (x, k x ) - S T &-f 9 Sy/ST (x, fe ± )) (Re(M!M*) cos <p* - Im^M^) si 



S L Afl /SL (x, kx) + S T &fl ls ^ (x, k x )) (lm{MiM%) cos <fi + Re(M L M*) sin V h q 



(80) 



A N D h/q ,(z,p x ) 
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Notice that this expression, at leading twist, is exact at all orders in k±/Q. We list here some properties of the 
polarized distribution functions which are useful in computing the asymmetries [2]: 

fq/S T (x,kx) + fq/-S T (x> k±) = 2f q/p (x,k ± ) 
f q /s T (z, kx) - f q /-s T (z, k±) = A/ 9 / St (a;, fej_) 
Af Sl /s T (x,kx) = -Af Sx/ _ ST (x,kx) 
Af Sy/ST (x,kx) - Af Sy/ _ ST (x,kx) = 2A-fl /ST (x,kx) 
Af Sy/ST (x, k ± ) + A/ SH/ _ Sr (x, kx) = -2 Ap sy/p (x, k±) 

Af Sz /s T (x,kx) = -Af Sz/ _s T (x,kx) 
Af Sl /s L (x, kx) = AJ Si /_ Sl (x, k ± ) (i = x, y, z) . 



(81) 



Let us now consider Eq. (80) in several particular cases. In the sequel, transverse and longitudinal always 
refer, both for the protons and the leptons, to their (different) directions of motion in the 7* — p cm. frame. 
Longitudinal states coincide with helicity states. 



A. Nucleon transverse single spin asymmetry, Aut 



Let us start with one of the most common SIDIS single spin asymmetries, Ag e s, with unpolarized leptons 
(U) and transversely polarized protons (T): 



(]6 a tpt^e'h.X _|_ ^6 a lp^^e'hX ~ J6 a i+p(S T )^l-'hX _|_ '^i+pJ^S^^hX 

For the numerator of Aut we have: 



dx B dQ 2 dzh d 2 Pr dips 



2vr ^ 16tt(x b s) 2 



d 2 kx — J 

Zh 



x \ -Af q/ST {x,kx){\M l \ 2 + \M 2 \ 2 )D h/q {z,px) 



A 



(82) 



(83) 



-f q Sy/ST (x,kx) (Re(M!M 2 *) cos^ - Im^M*) sin^)) 
Afl /ST (x,kx) (Re(M x M 2 *) sin^ + Im^hM*) «*#))] A N D h/qf (z,p ± ) 



The first term in Eq. (83) corresponds to the Sivers effect, whereas the second and the third terms correspond 
to the Collins effect, coupled to the transversity distributions. 
Similarly, for the denominator we find: 



da i P ^fhx + d(J e P ^i'hx 
dx B dQ 2 dzh d 2 Pr d4>s 



Ly 1 

!t ^ 16n{x B s) 2 



2tt 



d 2 kx — J 

Zh 



x iyf q ,p{x,kx) (|Mx| 2 + \M 2 \ 2 )D h/q (z,p ± ) 

-A.f q Sy/p (x,kx) (Re^M*) cos ^-Im(M!M 2 *) sin^ 1 ) A N D h/q1 



(84) 



Here, the first term corresponds to the usual unpolarized cross section (which survives in the collinear limit) 
whereas the second term is an effect obtained combining the Boer-Mulders distribution function, Af q , (x,kx), 



with the Collins fragmentation function, A D h / q t(z,px)- 

If we insert the exact relations for M\ and M2 - given in Eqs. (51) and (52) - and for cos ip q , s'mifi q - given 



in Eq. (C3) - into Eqs. (83) and (84), we obtain an exact expression for the Aut asymmetry. As already 
mentioned, the numerator is given by two different contributions, the Sivers and the Collins effect. Similarly, 
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the denominator, which is simply twice the unpolarized cross section for the £p — > t'hX process, receive most 
contribution from the first term, proportional to the unpolarized distribution and fragmentation functions, with 
a further contribution from a combination of the Boer-Mulders and Collins effects. 

Much simpler, and often quite accurate, expressions can be obtained at 0(k±/Q), neglecting higher order 
corrections. Using Eqs. ([57V(l60| and (p6|h(|32|) in Eqs. (l83l and (pil), one has: 



T ->£'hX _ ^Ip^^i'hX 

dx B dQ 2 dzh d 2 Px dips 

^£ X>* / <P*± {l A fv/s T M<f>± - <Ps) [1 + (l - y) 2 ] D h/q 

+ (1 - y) (Afl /ST + A-^ /Sr )(sin(^ + fa) - z h ^ sin(0 ± + fa)) A N D 

(1 " y) {Af ^ s - ~ A "^/5 T )( sin (^ + 2 ^ - M - MH± - 4> S ))A N D h/q 

(2 - y)^V^{Afl /ST + A-^ /ST )(sin(^ - fa + fa) -z h ^- sinfa)A N D h/q , 



p± ^ y q K~j. a/8T 1 - J Sy /s T . 



+(2-y)y/l^y"±Af, 



_ 2a 2 



^ Af q/ST ( sin fa - sin(2^ ± - fa)^D h/q 



sin(fa + (j>j_ - fa) - Zhjr sm(2^j_ - fa)) A N D h/q r 



{[l + (l-y) 2 ] sm(fa-fa)F^ 



,sin(0 fe -0 s ) 



(85) 



sm(fa + fa) F ^ Ms) + sm(3fa - fa) F$M>>-M 
sin 



and 



+ 2(1-2/) 
+ 2(2 -y)^l^y 

dx B dQ 2 dzh d 2 Pr dips 



An fa + S m(2fa - fa) F ^~^ } 



2q 



14 5Z e <? 
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/p 



J d 2 k ± |/, /p 1 + (1 - y) 2 - 4(2 - y) /T^y ^ 



COS 0j_ 



D 



h/q 



(1 - y) ( cos(fa + <j>±) 



k± 1 
z h -5- cos(20jj 



-2(2 - y) \j\-y — (cos fa - z h — cos fa_ 



2a 2 ( 

= _ | [l + (1 - y) 2 ] + 2(1 - y) cos 2^ 

where we have also exploited the definitions of the F structure functions, Eq 
) and ( 86 ) , can also be obtained directly from Eq. ( 79 1 . We recall 
) 1 = Pt — Zh k 1 and J = 1 . 



(86) 



2 ^+2(2-y)yr^cos^^^ 

structure functions, Eqs. (63|-(78). These last expressions, 
)m Eq. (79 1. We recall that, at 0(k±/Q), one has x = x B , 



VV AA^J. \_, VV^- 11UV\_. LllOVJ ^JYJJIWI Ul^U Lj A AVj ^IVjAA 

Eqs. (85) and (86), can also be obte 
z = Zhi P± = Pt — Zhk± and J = 1. 

The first term in Eq. (85 ) corresponds to the SIDIS Sivers asymmetry, which we analyzed in Refs. [i34. f40l - l42] 
for the extraction of the Sivers function, while the second term corresponds to the SIDIS Collins asymmetry, 
studied in Refs. H3] and used for the simultaneous extraction of the Collins and transversity functions. 



B. Nucleon longitudinal single spin asymmetry, Aul 

This asymmetry is defined for unpolarized leptons and a longitudinally polarized proton target: 

^> a lp^^l'hX _ fl6 a ip^^rt'hX ( j6 a e+p(S L )^e'hX _ a l+p(-S L )^l' hX 
^ UL = fffjlp^^i'hX + ( ]6 (T t P ^-+i'hX ~ dGpl+piSL^l'hX + dG a l+ v (-S L )^l'hX • ( 87 ) 
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We give explicit results, for this and the next asymmetries, only valid at 0(k±/Q). The denominator, as in the 
previous asymmetry, is twice the unpolarized cross section and is given in Eq. (|86[). For the numerator we have: 



da e+ P (s L )^'hx _ da e +P (-s L )^e'hx Aa z 



as can be easily checked from Eq. ( 79 1 



C. Nucleon longitudinal double spin asymmetry, All 

This asymmetry is defined by keeping fixed the longitudinal polarization of the lepton, while flipping the 
direction of the proton target longitudinal polarization: 

d 6 a e^ P ^^e'hx _ d 6 a e^p^^£'hx d e a e(s e )+ P (s L )^thx _ d e a i(s e )+ P (-s L )^e'hx 

^~ LL ~ d <o n t^p^^l>hX + d <o a £^p^^i>hX ~ d G a t(Se.)+p(S L )^t'hX _|_ d 6 a t(Se)+p{-S L )-^e'hX ' *- 89 ) 



The denominator is the same as given in Eq. (86), while for the numerator we have 

da t(S e )+p(S L )->e'hX _ d(J t(S t )+ P {-S L )^l'hX 



dx B dQ 2 dzh d 2 Pr dips 

^ {[1 - (1 - y) 2 ]F LL + 2y^T^j cos<f> h F™* h + 2(1 - y) sin2^F^ + 2(2 - y) ^/T^y sm cj> h F$l*») . 

(90) 

D. Lepton longitudinal double spin asymmetry, All 

This asymmetry is defined by keeping fixed the longitudinal polarization of the proton target, while flipping 
the lepton longitudinal polarization: 



d 6 a t-'p-*->t'hx _ d s a i^p^^e'hx d 6 a £(Si)+p(S L )^e'hX _ d 6 a l(-S e )+ P (S L )->l'hX 
All = ( jp (7 l-+p-*^i'hX _|_ d 6 a i^ P ->^i'hx = d 6 a e(s e )+ P (s L )^e>hx + d 6 a .e(-s e )+ P (S L )^e'hx 

For the numerator we have 

da e(s e )+p(s L )^e'hx _ da i(-s e )+p(s L )^i'hx 2a 2 



dx B dQ 2 dzh d 2 PT < 



{[1 - (1 -y) 2 ]F LL + 2y^y cos^F™ 8 ^} . 



(92) 



Notice that, in this case, the denominator differs from that given in Eqs. (86), as it acquires additional terms 
(generated by A/^ /s J: 

da e(s e )+p(s L )^£'hx + da i(-s e )+ P (s L )^i'hx 
dx B dQ 2 dzh d 2 PT d(j)s 

2 " 2 ( [l + (l-y) 2 } F uu+ 2(l-y)[cos2<p h F^+sm2<p h F^} 

+ 2(2 - y) ^~y [cos <j> h F^ *» + sin <f> h F*%+»] ) . (93) 



E. Nucleon longitudinal-transverse double spin asymmetry, Alt 

This asymmetry is defined by keeping fixed the longitudinal polarization of the lepton, while flipping the 
proton target transverse polarization: 

d 6 a e^p^£'hX _ j6 a £- t p i ^t'hX d e, a l(St)+p(ST)^'hX _ d 6 a t(S e )+p(-S T )^i' 'hX 
^ LT = d 6 a >t->pt->i>hX _|_ d 6 a l->pi—>e'hX ~ d 6 a l{S e )+p(S T )^t'hX _|_ ^6 cr £(5 f )+p(-S T )^<!'/iX ' ( 94 ) 



1G 



The denominator is given in Eq. ( 86 1 , while for the numerator we have 

da e{St)+p(s T )-yi'hx _ da e(s e )+p(-s T )^e'hx 
dx B dQ 2 dzh d 2 Px dips 



2a rr 

or{[i + d-»)- 
+ [l-(l-y) 2 
+2y^l - y cosfaF™ 
+2(1 -y) 



sm(fa - (ps)t UT 



cos(fa - s )Fr T s( ^ 0s) 



0S +cos(2 ( /», l - ( /» s )F^ s(2 ^- 0s) 

■sin(30h-0s) 



+2(2-y) v /l-y cos<?!> 



sin(^ + <fe)F™ ( ^ +fe) + sin(3^ - fa)F^ : 



, n j- 1 COs(0 h -0 S ) 



-■ os ^+co S (fa-fa)F^ 
sm <p s F UT 



7 iCua 1 

COS (ps-r LT 

psinips , _• In a 1 \ T-,sin(20 h -0 s )l 1 



(95) 



F. Lepton longitudinal-transverse double spin asymmetry Alt 



This asymmetry is defined by flipping the direction of the long: 
fixed the proton target transverse polarization: 



A - d&<7 



■ngitudinal polarization of the lepton, while keeping 

(96) 



l^pf^l'hX _ Mgl^pt^-l'hX c j6 a .e(Se)+p(S T )^t'hX _ ^,^.(-S t )+p(S T )^-£'hX 



For the numerator we have 

da e(s e )+ P (s T )->Chx _ da e(-s e )+p(s T )->i'hx 



dx B dQ 2 dzh d 2 Px < 
^ {[1 - (1 - yf] cos(fa - fa)F™^~^ + 2y^T—y 



cos 



faF™* s + co S (2fa - fa)F c L p 2<t>h " t ' s) ] } . 



The denominator differs from that given 
in the numerator of Aut- 

da l( Se )+p(S T )^i'hX + da 
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(97) 

Eq. (86), as it acquires several additional terms, which also appear 

£(-S e )+p(S T )^rt'hX 



2a 2 

Q A 



dx B dQ 2 dzh d 2 Px dfa 

sm(fa-fa)F^-^ ) ] 



|[l + (l-y) 2 ] [F uv 
+ 2(1 - y)[cos 2 fa F v °^ h + sm(fa + 

+ 2(2 - y)y/l=i [cos fa F^" + sin fa F^ s + sin(2^ - fa) F$F* h -+ s) ] 



fa) F v %^ s) + S in(3fa - fa) F** 3 **-*^ 



(98) 



G. Other asymmetries 



All the other single and double spin asymmetries are either zero or related to those already shown above. In 
particular, all the single spin asymmetries generated by the lepton polarization vanish: Ajjj — as Fjjj = 
to leading order in k±/Q and Atu — as we have no access to the transverse polarization of the lepton and 
therefore there are no terms proportional to either P^ or Py in Eqs. (jZJ) or (79). For the same reason we have 
A TT = A UT and A TL = A UL . 
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IV. PHENOMENOLOGY OF SPIN ASYMMETRIES 



To leading order in (fcj_/Q), all terms contributing to the polarized SIDIS cross section and to the spin asym- 
metries can be integrated analytically, provided we adopt a simple k± and p± dependence for the distribution 
and fragmentation functions. As usual, we assume the x and fcj_ dependences to be factorized and we assign the 
k± dependence a Gaussian distribution with one free parameter to fix the Gaussian width. For the unpolarized 
and helicity distribution functions and for the fragmentation function we simply use 

e -k 2 ± /(k 2 ± ) 

f q / P (x, k±) = f q/p (x) ^ (99) 

e -kl/(k 2 ± ) L 

Afl JSL ( x ,k ± ) = Afl /SL (x) (100) 



e -pi/<pi> 



D h /q{z,Px) = D h/q{z) _/„ 2 \ ' ( 101 ) 



where f q / p (x), A/* /g L ( x ) an< ^ ^h/ q ( z ) can De taken from the available fits of the world data. In general, we 
allow for different widths of the Gaussians for the different distributions, but take them to be constant and 
flavor independent. For the Sivers and Boer-Mulders functions, we assume a similar parametrization, with an 
extra multiplicative factor k± to give them the appropriate behavior in the small k± region (40] : 



— e xl c- 



M s 7r(fc 



2 



A/* ,„(*,*,) - A/! , (x)V2e 



k± -k 2 /m 2 e- fc i/< fc i> 



where the x-dependent functions Af q /g T (x) and A/' are not known, and should be determined phe- 

nomcnologically by fitting the available data on azimuthal asymmetries and moments; the k± dependent Gaus- 
sians have been assigned a reduced width to make sure they fulfill the appropriate positivity bounds: 



(fci) + M| 



(ki)s = Tcrrrik ( 104 ) 



{k2± )bm - jm+w ' (105) 

> -LI BM 

Similarly, for the distribution of longitudinally polarized quarks inside a transversely polarized proton, Af^ , 
and of transversely polarized quarks inside a longitudinally polarized proton, A/* , Sl , we set: 

= *z.,*rto^jt-^8r- (106) 

, &• , ,2 ;,,2 P~k±/(k ± ) 

, fci P~ k l/( k ±) TL 



with 



- MyrSt (108) 
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(k 2 ± ) M 



2 

TL 



( ^ >Ti {kl)+NP T 



(109) 



For the transversity distribution function, it is most convenient to parametrize the following combinations 



-ki/{kD T 



(k 



2 ^ 

_L/t 



(110) 



A/ s 9 x/ST (x,fc ± )-A-/^ /ST ( a; ,fc ± )) = 



2M 2 



hi T (x, k±) — h 1T (x) 



ek\ -k 2 x /M 2 e fc -L/< fc ^> 



AT 2 

TT 



p k 2 p (^±)q 



M 2 



AkV, 



(111) 



as these are the quantities which appear in the polarized cross section and in the spin asymmetries. Notice that 
for hi(x, k±) and h^ T (x, k±), as for each of the other TMDs, we introduce their own reduced Gaussian widths 



(k±} T 



(fci)M, 



TT 



^ TT (k 2 ± ) + NP 



Finally, for the Collins fragmentation function we choose 



(112) 



= A N D h/q Az)V2e 



P± e 



having defined 



(pi) M l 
(pD + m* 



(113) 



(114) 



Using the parametrizations in Eqs. ( 99||114 ) we can perform the k± integrations analytically in Eqs. ( 63|78 ), 
and re-express all the F structure functions in terms of the Gaussian parameters: 



F 



uu 



ycos 2<f> h 
UU 



pcos <p h 
r UU 



_p2 //p2 

e Ul/ P ( X B) D h/q(Zh) , p2N 



Pt 

Q 



M q Sy/p (x B ) A N D h/q ,(z h ) ^^/(^U Z h (fc 2 ) 2 M (pj) 

(ki)(pl) 



Af c 



P —Pt/( p t) 



A/^OO A N D h/q t(z h ) e^/i^BM 



M, 



t)bm 



b UL 



b UL 



b LU 



(kl)(p\) 

Mt /SL {x B ) A N D h/q r(z h ) e i-^/(^> M z^fcD^H) 



!(fci)BM(^-(^) B «)+^i) c ^ 



AL 



E 



A/*. /Si (x B ) A"B h/ ,t(%) e i-^/(^> Ti 



AL 



Jfcj)L(^i) 2 - 
(fci)H) 

at leading twist 



M h 



*(Pt) A 



T I TL 



l(ki) TL (P T ~(PT) TL )+(P 



Jc(Pt)tl 



(115) 
(116) 



(117) 



(118) 



(119) 
(120) 
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LL 



TpCOS lp 
f 1 



where 



v-^ » - - - - , - e p t/( p t) l 



< A fl/s L (*B)D h/q (z h ) 2 (121) 



KZ* h = -2^E ^/5,fa)^/«W e ~7<ft)\ Zh{klh (122) 

-^Z^ e ? jTf D h/ q (zh) — -r52y2 TVTV- ( 123 > 



'„2 \2 



pam&H-M p ST ^/St^b) , , e P t/(P 2 t) lt z^fcjj- 



'/ 



T/lt \ n X/ 

,N n . , *t*. \ „\li-Plnpl)„ /„2\2 



ptoWt-M P't y p2 fajr(^) e 3 ^'" jM^l^ fmi 



1 1 _ 2 A^fe) i/2-p-/(p- )s (fc^((pj)(P^) s +zgpj(fcj) s ) 

TlQ V e? D h/q( Z h) <fc2) (129) 

32 u± („ \ AN) 



F sin(20 h -0 S ) _ _ /k V e 2 A D h/qt {z h ) 



X 



3/2-P|/(P=) Tr ^D^Dc [^<fci) TT (P| - (P|) TT ) + 2K) C (P|) T 



_L/ 



75 g- 2. e « Ms Dh ^ Zh) ff <j?>3 ^iT (130) 



<p|> = (pD + ^^d 

(P 2 ), = (pD + ^^i)/ (I = S,L,LT) (131) 
(P|)j = (pi) c +z 2 (fci), 7 (J = T,BM,TL,TT). 



The unpolarizcd SIDIS cross section and all the asymmetries presented in Section III can now be rewritten 
in terms of the Gaussian-integrated P's, which depend on the TMDs. In order to single out information on 
a particular TMD from the measurements of the asymmetries, one has to disentangle the different azimuthal 
dependences. For example, the unpolarized cross section, see Eq. (86), includes the usual unpolarized collinear 
SIDIS cross section, the Cahn effect proportional to cosfih (studied in Ref. [34]), and a contribution generated 
by a combined Boer-Mulders and Collins effect, which appears in terms proportional to cos2</>h and coscj)h- 
Similarly, in the numerator of the Aut single spin asymmetry, Eq. ( |85[ ) , the Sivers and Collins effects are both 
simultaneously at work, together with other azimuthal modulations. To extract single effects, one introduces 
appropriate azimuthal moments of the asymmetries, defined as 

[da i(s e )+ P (s)^'hx _ da i(s t )+ P (-s)^hx ] w ^ 0s) 
h d<}>s [do*&)+p(S)-+t?hx + da e ^ + ^- s ^ thX ] 



Alt^ = 2 1 - f , (132) 
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where the function W{4>h,4>s) is an appropriate "weighting phase" which, upon integration, singles out one 
individual term of the asymmetry. For instance, to isolate the Sivers effect one can consider the sm((j>h — <fis) 
azimuthal moment of the Ajjt asymmetry: 

dfa dfa [dcT £ P r ^' hX - da e P^ e ' hX ] sin(fa - fa) 



A sin(0 h -0 S ) = 2 j (133) 

The W weight selects the Sivers term of the asymmetry in the numerator, while the integration over the az- 
imuthal angles <fis and <fih leaves only the first term of the unpolarized cross section, Eq. ( 86 ) , in the denominator: 

thus, this azimuthal moment is simply proportional to the ratio J F^ T / f Fjjjj. 

Furthermore, experimental data deliver these azimuthal moments as a function of one variable at a time, 
either x B , z^ or Py. Therefore, one has to integrate the numerator and denominator separately over all variables 
but one, in order to obtain the appropriate expression to be compared with the data. Clearly no simplification 
of common terms in the numerator and denominator can be made before the integrations have been performed 
(notice also that y is a function of both x B and Q 2 ). 

Let us consider, as an explicit example, the Sivers azimuthal moment A s ^^ h ^ s '(z/j), as function of Zh alone. 
Using the Gaussian-integrated expression of p^-^ h ~'t' s ^ f Eq. (f23) and integrating analytically over Pt we 
obtain 

f 

1 sin(0 h -0 s ) / _ ^ " ^ i 

i + (i - y) 2 

dx B d Q 2 Qi e'g f q / P {x B ) D h/q (z h 



dx B dQ 2 1 + V)2 £ e 2 q A N f q/ST (x B ) D h/q (z h ) 



A^-^\z h ) = A s — ,^ n _J^_ , (134) 

z 



where A s is a factor which only depends on Zh and on the free parameters which give the Gaussian widths for 
the distribution and fragmentation functions 

/ o „_ /J.2 \2 

(135) 




Notice the further dependence on zh hidden in {P T ) S) Eq. (131). 

Repeating similar procedures one can extract information on the other TMDs. The a zimuthal moment 
Af T iMs \ obtained using the weighting phase W(4>h, <f>s) — sin(^>/t + <fis) m Eq- (|132|) with unpolarized 
leptons, selects the Collins effect, coupled to the transversity distribution Ffz(x) = Atq(x) = hi(x). In this 
case, the azimuthal moment is sensitive to the ratio F^ T ^ h+ ^ s ' > / Fjjjj , and precisely: 



J dx B dQ 2 ^f^ ^2e 2 q h 1 (x B )A N D h/q ,(z h ) 



A s ^ Ms \z h ) = A c — « , (136) 

^2 1 + I 1 ~ V) \ " 2 

Q4 e 8 • 



dx B dQ 2^e q f q/p {x B )D h/q (z h ) 

i 



with 



A °"4M^V^W ' (13?) 

One can further exploit the Ajjt asymmetry, to isolate and measure the transverse distribution function 
F7j(x) = h^ T (x), by weighting the single spin asymmetry numerator with the phase W(4>h, 4>s) = sin(30/j — 4>s), 
obtaining: 



dx B dQ 2 2{ ^ J 1 £ el hi T {x B )A N D h/q ,{z h ) 

4 sin(30 h -0 s ) /„ \ _ a "_ 1 /i qo^i 

J dx B dQ 2 — ^e 2 q f q/p {x B )D h/q {z h ) 
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where 



3ezg 2 en (kj)^ {pj)l 

8M* T M h (p2) TT V (P*) TT (kl) (pi) 



A - Z ' e " \_-L/_TT \T±JC / 13 Q\ 

tt a n.r2 \/f. /d2\ \ /p2\ ll„2\ /„2 \ 



One can write similar expressions for all other asymmetries, which we do not report here. From A s ™jj h 
and A s ™ 2< ^ h one can obtain information on Af Sx /s L , while A^^ s , A c ^^ h and A™y ( ' 2 ^ h ^ depend on 
Af Sz /s T - A^f s and A s ^, 2 ^ h are more complicated to analyze as they receive contributions from the Sivers 
distribution function (both of them) and, in addition, from the transversity distribution h\{x) (A^^ s ) and from 

(.1 I Asin(2tph-<t>s)\ 

"it \ a ut )■ 

Let us consider in more details the unpolarized cross section, to which, remarkably, a similar "weighting'' 
procedure can be applied. In fact, one can introduce the average value of W(<ph) with an expression similar to 



Eq. ( 132 ) in which the unpolarized cross section appears in the numerator as well as in the denominator 



dd> h d^ s [d*^' hx + da^' hx ] W{4> h ) 

(Wtfh)) = J - f ■ (140) 

d</> h d^ s [d«A Wwc + da*P l ^ hx ] 



For instance, weighting the unpolarized cross section with W(cj>h) = cos20^ one can gain direct access to the 
Boer-Mulders function, coupled to the Collins function (on which independent information can be obtained): 



J dx B dQ 2 V^^elAf« y/p (x B )A N D h/qf (z h ) 



(cos2 0/l ) = A BM — 2l + {l l y) 2 ~ , 



dx B dQ — A 2^eqf q / P {x B )D h/q {z h ) 



with 



ezu (kj) 2 BM (pDI 

BM M BM M h (P*) BM {k\) (pi) ' [ } 



Analogously, using W(<f)h) = cosfih, one has 



J dx B dQ 2 (2 E e l [ A -p M* B )Dh/M) + B BM Af q Sy/p K)^ N D h/q ^z h ) 
= " \ 2 l + (l-,) 2 ^ 2 ,U:!1 



dx B dQ — A 2^e q f q/p (x B )D h/q (z h ) 

g 



with 



^unp — *h r\ \ IT>2\ ' d bm — r,^l»/r Tilt ll-2\ /„2 \ , „n,3/2 



unp Q V^l> ' BM 2 Q M BM Mh (kl) (pi) (P 2 } 3/2 
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V. CONCLUSIONS AND FURTHER REMARKS 

The study of the 3-dimensional structure of protons and neutrons is one of the central issues in hadron 
physics, with many dedicated experiments, either running (COMPASS at CERN, CLASS at JLab, STAR 
and PHENIX at RHIC), approved (JLab upgrade) or being planned (ENC/EIC Colliders). The transverse 
momentum dependent partonic distribution and fragmentation functions, together with the generalized parton 
distributions, play a crucial role in gathering and interpreting information towards a true 3-dimensional imaging 
of the nucleons. TMDs can be accessed in several experiments, but the main source of information is semi- 
inclusive deep inelastic scattering of leptons off polarized nucleons. The theoretical framework in which the 
experimental information is analyzed is the QCD factorization scheme. 

We have used here an intuitive approach to TMD factorization in SIDIS and shown that one can re-derive, at 
leading order, the most general expression of the polarized cross section, obtained within the QCD factorization 
scheme by other authors [2 Q20 [22] ■ All azimuthal dependences are precisely generated by the properties of the 
helicity amplitudes, which we use to describe the factorized steps of the process: the partonic distributions, the 
elementary interaction and the quark fragmentation. 

We have obtained explicit expressions for all the SIDIS spin asymmetries and the cross section azimuthal 
dependences which allow to extract information on the TMDs. Indeed, some of them have already been used 
to study the Sivers [35] [35], the Cahn [IS] [33] and the Collins [3J effects. Simplified expressions, based on a 
Gaussian k± and p± dependence of the distribution and fragmentation functions, recently supported by data 
[4 7) , have been given; they might be useful for fast and simple analyses of the experimental data. 

We wonder, at this stage, whether the same approach can be used for other processes. It works, with the 
same validity as for SIDIS, for Drell-Yan processes (D-Y) [33], where our helicity amplitudes for the different 
factorized steps reproduce the most general azimuthal structure of the cross section as obtained in the TMD 
factorization [TS]. As commented in the Introduction, both in SIDIS and D-Y the presence of two different 
natural scales, a small and a large one, is crucial for the validity of the QCD TMD factorization. 

Our approach was actually first introduced for processes with a single large scale, like pp — > irX, with large 
Pt pions [2]. These are the processes for which the largest single spin asymmetries have been observed and 
might be generated by TMDs 48 50 . However, TMD factorization has not been proven in these cases. Despite 
that, an extension of the intuitive approach used for SIDIS - and shown to be perfectly equivalent to the QCD 
TMD factorization scheme - is natural. That was the guiding idea in Ref. [2]; each proton "emits" a parton, 
the two partons interacts and one of the final parton fragments into the observed hadron. All intrinsic motions 
are taken into account and phases appear in the helicity amplitudes. The difference with SIDIS processes is 
that, in this case, the measured large Pt of the final hadron is generated by the hard elementary scattering, 
and all intrinsic motions are integrated over. As a consequence, the phase integrations strongly suppress the 
relevance of most TMDs, with the exception of the Sivers and Collins effects [51] [52], which combine into the 
observed asymmetry, and cannot be separated unless one could resolve the internal structure of the final jet 

m 

A global simultaneous phenomenological analysis of single spin asymmetries in SIDIS and pp interactions 
is, at the moment, rather difficult. Apart from the validity of the factorization scheme in both cases, another 
important open point is the universality of the Sivers functions; it is not clear whether or not they should be 
the same in the two processes or should be corrected by some gauge color factors [30]. In any case it is worth 
trying to explore the possibility to have a unique description of SSAs in different processes, based on TMDs; 
work in this direction is in progress and will be presented elsewhere. 
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Appendix A: Helicity Amplitudes 

We show the explicit computation of the helicity amplitudes M\ 3 x i - t \ 1 x 2 f° r the non-planar process i(k%, Ai) + 
<z(&2, A 2 ) — » H'{k 3 , A 3 ) + <7'(fc4, A4), in the 7* — p cm. frame of Fig. 1. We exploit the spinor helicity technique, 
adopting the conventions of Ref. [35] ■ At LO in QED, when neglecting all masses, there are two independent 
helicity amplitudes: 

= (q' + \rk + ) (e' + \l,\t + ) = (4 + | 7 l2+) (3 + | 7ai |1+> (Al) 

= (r\r\q-) = (4-| 7 ' 1 |2-) (3+| 7ai |1 + > , (A2) 



which can be written as 



where 



M ++ , ++ = 2 ^-[43] (12) (A3) 
M+-.+- = 2 ^ [23] (14) , (A4) 



u_(fc i )tt+(* i ) = (ij) = -fol* = ^/kjkje-^-^/ 2 - JkTkfe^-^/ 2 (A5) 

u+(fc*)M%) = [y] = -<y'r. (A6) 



with fc± = fc° ± k 3 . 

In the 7* — p cm. frame we have (see Ref. [53] for details): 

fci = £7(1, sin^, 0, cos 6) 

If Q 2 Q 2 

k 



k 3 = ki - q 



k 4 = k 2 + q 
01,3 = , 02,4 

where, neglecting the proton mass: 



1 / / k 2 
x = -x„ I + 4/I+4 - 



E 



2 ~B ■ y- ■ - Q2 

s-Q 2 y/s l-x B y 



2W 2 y/y(l-x B ) 

2 



x B ys 

P 



W=^y{l-x B )s 



cose = 1 + ( " - 2 »' T - si„o = VjgZKBZg 

1 - yx B 1 - yx B 

These relations allow us to express all the kf components in terms of x B and y [3 



k+ = E{l + cos6) = y/s< '' 
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k 2 



k^ 



E(l - cost 
E(l + cos( 



; x B {l- y) 

VVU- - X B, 



Q J — -^n 



E(l -cos (9) - — = Js- r 
1 ' W ^y^xj) 



2xPq Xy/S 

2xP = x\fs 



y 

y 

1 — x„ 



2xP, 







fc7 = 2xP - 



Q 2 



1 - x D 



ll 



V 



t>3 = , 02 = 04 = 0_L • 



From Eqs. (|A3])-(|A6) we get: 



^1 ^2 V ^1 ^2 e 



^3 ^4 \/ ^3 ^4 ^ ~* 



M. 



+-;+- 



= 2- 



^1^4 V^'l ^4^ e ^ J 



Exploiting Eqs. (A9| we can finally compute the amplitudes as function of y, Q and A;j_ 



M++-++ = e q e 



1 



l+4-| | < 



-i<t>± 



y 



1- 4/1 



fc2 



y Q 



(A9) 



(A10) 
(All) 

(A12) 



M+-.+- = e q e 2 



1 + 4 1 



fc2 

4 7# 



1 - 4 1 



fc2 

4 <? 



yi - 2/ fcj- 

2/ Q 



(A13) 



Appendix B: Helicity formalism and helicity transformations 

All our analytical and numerical computations of the SIDIS cross section, Eq. Q, are performed in the 7* — p 
center of mass frame (cm.), with the kinematics represented in Fig. [I] However, in our helicity formalism all 



components of the polarization vectors (like in Eqs. (17) and (18l) and of the transverse momenta which enter 



the definition of the TMDs, refer to the appropriate helicity frame of the corresponding particle. Then, in order 
to perform our calculations, we have to express the helicity frame variables in terms of the cm. ones, which 
requires some care. 

For the proton, which moves along — Z cm , the helicity frame (X p ,Y p , Z p ), as reached from the 7* — p cm. 
frame, is given by (as discussed in Appendix D of Ref. [2]): 



so that 



X p X cm Y p Y cm Zp Z cm , (1^1) 



k ± = cos ip ± X p + sin Lp ± Y p = cos (p ± X cm + sin cf) ± Y cm = cos ip. x X cm - sin (p ± Y cm 

fe 2 = k x - [x B P Q - ^p^j Zcm (B2) 

St = cos ip s X p + sin ip s Y p = cos fis X cm + sin 4> s Y cm = cos ip s X cm - sin tp s Y cm , 
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which implies <p± t s — 2n ~ 4>±.s- As long as there is no ambiguity we use if for angles defined in the helicity 
frames and <j> for angles defined in the cm. frame, following the notations of Fig. 1. 

It is less straightforward to deal with the quark polarization vector, P q = (P q ,P q ,P q ), which describes 
intrinsic properties of the proton constituents, and is defined in the quark helicity frame. In order to keep the 
same definitions, through the helicity formalism, of the polarized TMDs as in Ref. [3J, we have to define P q in 



the quark helicity frame as reached from the proton helicity frame, 
frame are then given by [2j [35] : 

Zq = k 2 

y q = Z p xkx = -Z cm x k± 



The axes x q ,y q ,z q of the quark helicity 



y q x z q = (Z p x k±) x k 2 = -{Z cm x k±_) x k 2 



(B3) 
(B4) 
(B5) 



Notice that the quark helicity frame as reached from the cm. frame (Z cm ) is different from the quark helicity 
frame as reached from its parent proton helicity frame (Z p ); although the z q axes obviously coincide, x q and y q 
have opposite signs, Eqs. (B5) and (B4). Therefore, when referring to the kinematical configuration of Fig. 1, 



which we use throughout the paper, we have to take the x and y component of the quark polarization vector, 
P£ and P£, with opposite signs with respect to those obtained from Eq. (15); this has been done in Eqs. (24) 
and p5l. 



Appendix C: Analysis of the fragmentation process 

Let us now focus on the azimuthal angle (p q involved in the fragmentation process. This is the azimuthal 
angle of the momentum P^ of the final hadron around the direction of the fragmenting quark q, as defined 
in the quark q helicity frame, see Fig. [2j Notice that the fragmenting quark, in the 7* — p cm. frame, has a 
longitudinal component along the positive Z crn axis. Its helicity frame, as reached from the 7* — p cm. frame, 
is given by Ref. [2J: 

Z = &4 

y = Z cm x k ± (CI) 
x = y x z , 

where fcj_ is the unit transverse component - with respect to the Z cm direction - of the outgoing quark, k<±. 

In the quark helicity frame, ip q coincides with the azimuthal angle which identifies the hadron transverse 
momentum , therefore 

cos if q = p ± ■ x 

sm<p*=p x -y. (C2) 




FIG. 2: Kinematics of the fragmentation process. 



2G 



By using the SIDIS kinematics as reported in Ref. [51] , one 

^z 



\ finds 



[P T kf cos{<j> h - <j>±) - k x ] 



sin ip q 



o\ = — S m(<t> h - £l) , (C3) 
P± 

where the superscript Z refers to the 7* — p cm. frame, where one measures Pt = (Pt cos <j>h, Pt sin <ph, P/f), 
and 




as derived in Ref. |34| . 

At (D(k±/Q) one simply has 



COS(/3 g = 



sin^ = 



Pt 
Px 
Pt 
P± 



COs((j>h - 4>±) - Zh^T 



sin((f) h - <j> ± ) , 



having neglected terms 0(k\jW 2 ) and 0{P T /W 2 ). 



Appendix D: Tensorial Analysis 



(C4) 



(C5) 



Eqs. ( 63 )-( 78 ) are obtained using a simple euclidean tensorial analysis, as outlined in what follows. In 



general, the tensorial structure of each of the P's functions defined in Eqs. ( 63 1-( 78 ) can be reduced to a linear 
combination of the convolutions 



T l = Jd 2 k ± Af(x,k x )kiAD(z )Px ) 
T 13 = J d 2 k ± Af(x,k ± )k i L k j L AD(z,p±) 
T jl = J cPk ± Af(x,kx)k ± k{k l ± AD(z, Pl _) 



(Dl) 
(D2) 
(D3) 



where we have denoted by Af (AD) any distribution (fragmentation) function appearing in the definition of 
the particular F function one is considering, while the fc^_, i = X, Y (X and Y refer to the 7* — p cm. frame, we 
have dropped the cm subscript) are the components of the k± transverse momentum vector, k^[ = fcj_ cos </>j_ , 
fcj_ = k± sin <f)± . One should bear in mind that p± is not an independent quantity, as it can be expressed in 
terms of k± and Pt- Notice that T\ T %] and T y are symmetric, rank 1, 2, 3 euclidean tensors respectively. 
Once the integration over d 2 k± is performed, the T l , T lJ and T rjl can only depend on the observable quantities 
Pt and (f>h, i-e. the measured modulus and azimuthal phase of the final observed hadron transverse momentum 
Pt- Therefore, in a completely general way, it must be 



T l = P T Si{P T ) 

T ij = P T P T S 2 (P T ) + S ij S 3 (P T ) 

T ijl = pi pJ pk 54 (p T ) + (pi gjl + pj gil + pl T § ij^ S - 5 (p T ) 



(D4) 
(D5) 
(D6) 



where the Pt components (Pj-f = P T cos(j)h, Pf = Pt sin <f>h) give the proper tensorial structure, while S1—S5 
are five scalar functions which can only depend on Pt (modulus), and can easily be determined by contracting 
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Eqs. (Dl )-(D3 1 with some symmetric tensorial structures (P T , , etc., as appropriate) to obtain simple scalar 
relations. Finally, one finds 

Si(P T ) = jrj d2k ±- ' P t) A/(x, k x ) AD(z,px) (D7) 

S 2(Pt) = ^2 J d2k ± [2( fc ^ ' P t) 2 ~ fci] A/(z, fc x ) AD(z,p ± ) (D8) 

S 3 (Pt) = / d 2 fe± [*i - (fc± • Pt) 2 } Af(x, AD(z,p ± ) (D9) 

•^(Pr) = ^3 / ^ [4( fc ^ ' P t) 3 - 3fcl(fe ± • P T )} Af(x, kx) AD(z,p ± ) (D10) 

S 5 (P T ) = ~ / d 2 fc ± [fci(fe x ■ Pt) - (fcx ' Pt) 3 ] A/fo fcj.) AD(*,;p x ) . (Dll) 
As a consequence, we have 

J d 2 kx cos cj)x Af AD = cos 4> h J d 2 kx (kx ■ P T ) Af AD (D12) 

J d 2 kx sin <j>x Af AD = sin <f> h J d 2 kx {kx ■ P T ) Af AD (D13) 

J d 2 kx cos 2 fix Af AD = ^J d 2 kx [l + cos20 ft [2(fc ± • P T ) 2 - 1]| Af AD (D14) 

y d 2 fc_L sin 2 4>x Af AD = ~J d 2 kx {l - cos 2<j> h [2(kx ■ P T f + 1}} Af AD (D15) 

/ A ± cos Si „ ^ A/ AD . cos , h si n / «V [2(t ± - P T >* - 1] Af AD (D16) 
y d 2 fc ± cos 3 <j)x Af AD = cos 3 t/> h y d 2 fc± [4(fe± • P T f - 3(ftx • P T )\ Af AD 

+ 3 cos 4> h J d 2 kx [{kx ■ Pt) - {kx • Pt) 3 ] A/ AD (D17) 
y d 2 kx sin 3 ^ A/ AD = sin 3 <j> h J d 2 kx [A(kx • P T f ~ 3(fc_L • P T )} Af AD 

+ 3 sin cj> h y d 2 fe ± [(fe_L • P T ) - (kx • Pt) 3 ] Af AD (D18) 
y d 2 fc ± cos 2 ± sin cj)x Af AD = cos 2 <j> h sin h y d 2 fc ± [4(fc_L • P T ) 3 - 3(fe x • Pt)] A/ AD 

+ sin <j> h J d 2 k ± [(kx ■ Pt) - (kx ■ Pt) 3 ] A/ AD (D19) 
y d 2 k ± cos sin 2 Af AD = cos h sin 2 0,, y d 2 fc ± [4(fc ± • P T ) 3 - 3(fe x ■ P T )] A/ AD 

+ cos 4> h y d 2 fc ± [(fe_L • P T ) - (kx ■ Pt) 3 } Af AD . (D20) 
From these equations one can easily reconstruct 

y d 2 kx cos 2cf)x Af AD = cos 2<f> h J d 2 kx [2(kx • Pt) 2 - 1] A/ AD (D21) 

y d 2 fc ± sin20_L A/ AD = sin 2^ J d 2 kx [2(kx ■ P T f - 1] Af AD (D22) 

y rf 2 fc ± cos 30j_ Af AD = cos 3<fo J d 2 k ± [4(k ± ■ P T f ~ 3(fc_L • Pt)] A/ AD (D23) 

y d 2 fc ± sin 30 ± A/ AD = sin 3<j) h J d 2 k ± [4(k x • P T ) 3 - 3(fe x ■ Pt)] A/ AD . (D24) 



All of these terms are easily recognizable in Eqs. (63l-(78) 
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Appendix E: Integration by rotation in the hadronic plane 



Eqs. (63)- (78 1 can also be obtained in a simple way looking to a slightly different reference frame. Let us 
define the production plane as the plane containing the virtual photon 7*, the proton momentum and the 
produced hadron h. We can define a new 7* — p cm. frame where the X' — Z' plane is the production plane. 
This new frame is rotated by an angle <f>h with respect to the cm. frame (X , Y , Z) depicted in Fig. 1 (we drop 
for simplicity the subscript cm): 



X = X cos <ph — Y sin iph 
Y = X sin tfih + Y cos <ph 



(El) 
(E2) 



Notice that X = Pt = h. Any integration in Eqs. ( 63 )-( 78 ) , at fixed values of the external variables, can be 
recast as the sum of one or more contributions of this kind: 



/d»*X*XOOB^/(*X,Px) 

J d 2 k ± k\ cos2<j>xf(kx,p±) 
J d 2 kx kx cos30j_ f(kj_,pj_) 



J d 2 kj_ kj_ sin <p±f(kx,Px) 

Jd^xkl^xfikx^) 

J d 2 kx k x sin30_L f(kx,P±) 



where 



pi =P4 + ziki-2z h (kx-p T )- 



(E3) 
(E4) 
(E5) 

(E6) 



Let us consider, for instance, Eq. (E3|; using Eq. (El), we have 



J d 2 k A 



xkxCOS(j)xf{kx,P±) = / d-kxkx f(k±,p±) = I d z kx {kx ■ X) f(k ± ,p± 



= d 2 k. 



(kx ■ X )cos</> ft - (fcj_ • Y )$m4>h /(fcj_,fcj_ ■ X ) 



cos(j) h / d 2 kx(kx- P T ) f(k±,p±) , 



(E7) 
(E8) 



where in the step (E7) we have underlined that / is a function of (kx ■ Pt) = (kx ■ X ) by means of Eq. (E6). 
With similar arguments we have, for all integrals of the kind (E3)-(E5): 

(E9) 
(E10) 
(Ell) 
(E12) 
(E13) 
(E14) 



d 2 kx kx cos0_l =>■ cos 0ft / d 2 kx (kx • Pt) 

d 2 kxkx^x^S^ h Jd 2 kx(kx-PT) 

d 2 k ± k\ cos 2(j>x cos 2cj) h J d 2 k ± [2(k ± ■ Pt) 2 - kj_] 
d 2 kx k 2 x sin20 ± =^> sin 2cf> h / d 2 kx [2(k x ■ Pt) 2 - fcjj 



d 2 k ± ki cos30_l cos3<t) h I d 2 k ± (kx • Pt) [4(fc_L • Pt) 2 - 3k 2 ± ] 
d 2 k ± k ± sin30_L => sin 3cf> h / d 2 kx (kx ■ P T ) [4(fex ■ Pt) 2 - 3k x ] , 



which coincide with Eqs. (D12), (D13) and (D21)-(D24). 
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